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Description 

Field of the Invention 

This invention relates generally to monitoring and 
analysis of the electrical activity of. for example, a 
patient's heart and, more particularly, to such monitor- 
ing and analysis of the spatial distribution of the electri- 
cal activity, based upon, for example, information 
collected from a limited number of electrodes. 

Background of the Invention 

Many physiological activities have electrical signals 
associated therewith. For example, the operation of the 
heart is regulated by electrical signals produced by the 
hearts sinoatrial (SA) node. Normally, the SA node pro- 
duces an electrical impulse roughly 60 times a minute. 

Each impulse generated by the SA node initially 
spreads across the atria of the heart. After reaching the 
atrioventricular (AV) node of the heart the Impulse then 
travels downward to the ventricles. The electrical 
impulse depolarizes the muscle fibers as it spreads, 
with atrial and ventricular contractions occurring after 
the impulses pass. After contracting, the muscle cells 
are unable to react again for a short interval of time 
known as the refractory period. Eventually, the muscle 
fibers repoJarize and return to their resting state. 

As will be appreciated, the electrical activity initi- 
ated by the SA node is, thus, representative of the oper- 
ation of a patient's heart lb allow the heart's operation 
to be analyzed, a variety of techniques have been devel- 
oped for collecting and interpreting data concerning the 
electrical activity of the heart. Perhaps the most basic of 
these approaches is the three-lead electrocardiogram 
(ECG). 

A three-lead ECG system typically employs a mon- 
itor and four limb" electrodes, attached to the patient to 
collectively monitor three voltages or leads." Specifi- 
cally, a left arm electrode (L) is placed on the patient s 
left arm. Similarly, a right arm electrode (R) Is placed on 
the patient's right arm. A left leg electrode (F) is 
attached to the patient's left leg and a right leg or 
"ground" electrode (Q) is attached to the patient's right 
leg. 

As an electrical impulse spreads across the heart 
the monitor repetitively measures the voltages at elec- 
trodes L. R, and F, relative to the ground electrode G. 
These voltages are designated vL, vR, and vF, respec- 
tively. 

The first of the leads 0) evaluated by the monitor is 
equal to the difference between the voltages at elec- 
trodes L and R ( l=vL-vR ). Similarly, the second lead (II) 
monitored with the electrodes is equal to the difference 
between the voltages at electrodes F and R ( H=vF-vR ). 
Finally, the third lead (III) evaluated by the monitor is 
equal to the difference between the voltages at elec- 
trodes F and L ( IH=vF-vL ). 

The waveforms produced by plotting each of the 



leads over an interval of time corresponding to one car- 
diac cycle are conventionally divided into a number of 
segments. The portion of each waveform representing 
atrial depolarization is referred to as the "P" segment or 
5 wave. Depolarization of the ventricular muscle fibers is 
then represented by "Q," M R," and "S" segments of the 
waveform. Finally, the portion of the waveform repre- 
senting repolarization of the ventricular muscle fibers is 
known as the T segment or wave. 

10 To simpGfy the analysis of the three leads I, II, and 
III of data, the heart is generally treated as a dipolar 
source of electrical activity, positioned at the intersec- 
tion of transverse, superior, and inferior axes of the 
patient's body. Each of the three leads Is regarded as 

is being equally sensitive to all sections of the heart Lead 
I is, however, interpreted as representing electrical 
activity on the transverse axis, while leads II and 111 are 
interpreted as representing activity on the superior and 
inferior axes. 

20 As will be appreciated, the leads are influenced by 
a variety of factors. For example, the size of the patient's 
heart, as weD as its location and orientation relative to 
the patient's trunk, will affect the magnitudes of the 
leads. The presence of heart disease will also typically 

25 impact the waveforms. Further, the size, fat content and 
air content of the patient's trunk as well as the exist- 
ence of any physical deformities or diseases, may influ- 
ence the waveforms. Finally, the relative positioning of 
the electrodes on the patient may alter the waveforms 
sensed. 

To properly interpret the three leads I, II, and III of 
data, a physician must be trained to recognize and iden- 
tify the influence of these various factors. The physician 
must also understand the dipolar cardiac model and the 

as resultant relationship between the three leads and the 
cardiac activity of interest. 

In that regard, a physician's analysis of three-lead 
ECG data typically begins with the selection of the par- 
ticular lead or leads to be reviewed. In that regard, 

40 because lead I reflects cardiac activity on the transverse 
body axis, the physician would most Gkely review lead I 
to detect lateral infarction or coronary vessel involve- 
ment. Similarly, because leads II and III reflect cardiac 
activity on the superior-inferior axes of the patient's 

45 body, they are reviewed to detect, tor example, inferior 
disease. 

Having selected a lead for analysis, the physician 
typically makes an Initial assessment based upon the 
overall appearance of the waveform. The physician 

so must also recognize the significance of particular fea- 
tures of the waveform. In that regard, a number of wave- 
form features have been identified for use in diagnosing 
various heart conditions. 

Tne area of the "QRST segment of the waveform 

55 has been recognized as indicating disparity of recovery 
properties. As will be appreciated, the waveform's area 
can be computed by Integrating the waveform over the 
Interval of time during which the ventricular muscle fib- 
ers are depolarized and repolarized. A change in QRST 
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area may indicate the presence of disease. 

Another waveform feature used to diagnose the 
heart's condition is an analysis of the "ST" segment of 
the waveform. In that regard, the displacement of the ST 
segment is generally indicative of acute injury. For s 
example, an elevated ST segment associated with Q 
waves indicates an acute or recent infarct. The analysis 
of the ST segment is usually based upon data collected 
from more than the three leads describes above. 

A related form of analysis is referred to as "late- 10 
potential" analysis. Late-potential analysis involves the 
expansion of the time-scale of the QRS segment of the 
waveform and a detailed review of the expanded seg- 
ment for features not apparent in a conventional display 
of the waveform. More particularly, the expanded seg- 75 
ment is reviewed for low level activity that persists into 
the ST segment. 

While the three-lead system described above is 
generally useful in evaluating the condition of the 
patient's heart, it does have certain limitations. For 20 
example, because the three-lead analysis is based 
upon a dipolar view of the heart, the leads have little 
selective sensitivity to the operation of particular regions 
of the heart. Depending upon the nature of the cardiac 
event to be detected, this lack of sensitivity may prevent 25 
the three-lead data from being diagnostically useful. 

One effort to increase the regional sensitivity of 
ECGs led to the development of "12-lead" ECGs. A 12- 
lead ECG system includes the limb electrodes of a con- 
ventional three-lead ECG system. In addition, a 12-lead 30 
system includes six "precordial" electrodes, positioned 
on the patient's chest. The six precordial electrodes 
were originally added to allow the electrical activity of 
individual regions of the heart to be more directly 
sensed. In practice, however, the dipolar model of the 35 
heart has been retained and the 1 2 leads interpreted to 
evaluate electrical activity along the three axes. 

Although the 12-lead ECG configuration has been 
described as allowing a limited form of "mapping" to be 
performed, for the purposes of the ensuing discussion it ao 
will be more precisely referred to as a technique for 
introducing selective sensitivity into the leads of an ECG 
system. In contrast, as will be described in greater detail 
below, mapping is more appropriately used to describe 
approaches involving the evaluation or presentation of 45 
spatial distributions of data over the patient's chest. 

Reviewing the construction and operation of a 12- 
lead ECG system in greater detail, as noted above, six 
precordial electrodes are added to the four limb elec- 
trodes of a three-lead ECG system. In that regard, a first so 
precordial electrode (p1) is attached to the front of the 
patient's chest, to the right of the sternum in the fourth 
rib interspace. A second precordial electrode (p2) is 
also attached to the front of the patient's chest, left of 
the fourth rib interspace. A fourth precordial electrode 55 
(p4) is positioned at the left midclavicular line in the fifth 
rib interspace. The third precordial electrode (p3) is 
positioned midway between electrodes p2 and p4. 
Finally, the fifth precordial electrode (p5) is positioned at 



the midaxillary line in the fifth rib interspace and the 
sixth precordial electrode (p6) is positioned midway 
between electrodes p4 and p5. 

The 12 leads or waveforms evaluated with the four 
limb electrodes and six precordial electrodes are as fol- 
lows. In addition to the three leads I, II, and III already 
described, three leads aVR, aVL, and aVF are formed 
from the voltages sensed at the limb electrodes. In that 
regard, lead aVR=vR-(vL+vF)/2 ; lead aVL=vL- 
(vR+vF)/2 ; and lead aVF=vF-(vL+vR)/2 . 

The six precordial electrodes p1, p2, p3, p4, p5, 
and p6 are used to sense voltages vl , v2, v3, v4, v5, 
and v6, respectively, referenced to the ground electrode 
G. The remaining six leads are designated V1, V2, V3, 
V4, V5, and V6. In that regard, lead 
V1=v1-(vR+vL+vF)/3 ; lead V2=v2-(vR+vL+vF)/3 ; lead 
V3=v3-(vR+vL+vF)/3 ; lead V4=v4-(vR+vL+vF)/3 ; lead 
V5=v5-(vR+vL+vF)/3 ; and lead V6=v6-(vR+vL+vF)/3 . 

As will be appreciated, the 12 waveforms produced 
by such a system each have the various PQ RST seg- 
ments described above. In analyzing the 12-lead data, 
the physician may consider the overall appearance of 
the waveforms, as well as the various waveform fea- 
tures discussed above in connection with the three-lead 
system. In addition, because the 12-lead ECG contains 
considerably more information than the conventional 
three-lead ECG, the physician may choose to review 
only a few of the numerous leads, with the particular 
leads selected depending upon the cardiac features of 
interest. 

In that regard, the dipolar model of the heart is 
almost universally relied upon in the analysis of 12-lead 
ECGs. Leads I, aVL, aVR, V5, and V6 reflect cardiac 
activity on the transverse body axis and are monitored 
to evaluate lateral disease. Leads II, III, and vF are 
indicative of electrical activity on the superior-inferior 
axes and are reviewed to evaluate inferior disease. 
Finally, leads VI, V2, V3, and V4 represent activity on 
the anteroposterior axis and are monitored to determine 
anterior and posterior disease. 

One alternative use of the conventional 12-lead 
ECG system described above is the performance of a 
stress test. The system is virtually the same with the 
exception that electrodes L and R are moved from the 
patient's arms to the trunk. As a result, leads I, II, III, 
aVR, aVL, and aVF are attenuated. Cardiac stress is 
introduced by subjecting the patient to a regimen of vig- 
orous exercise and the resultant leads are analyzed in 
conventional fashion. 

As noted previously, although the 12-lead ECG 
originally represented an effort to adopt a regionally dif- 
ferentiated model of cardiac activity in place of the dipo- 
lar model, both the three-lead and 12-lead ECG are 
conventionally analyzed as representing the activity of a 
dipolar cardiac generator. Serious efforts to provide 
selective sensitivity to the operation of particular regions 
of the heart were initiated by the introduction of body 
surface mapping (BSM) systems. Such systems evalu- 
ate or present spatial distributions of data over the 
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patient's chest. 

Conventional body surface mapping systems 
employ a relatively large number of electrodes in com- 
parison to the three-lead and 12-lead ECG systems 
described above. For example, much of the work done s 
in the area of body surface mapping has been with sys- 
tems employing 192 electrodes. As will be appreciated, 
the use of such a large number of electrodes allows the 
spatial distribution of voltage across the patient's chest 
to be relatively closely evaluated. to 

Unlike the ECGs discussed above, in which each 
lead is displayed as a time-dependent waveform, the 
voltages sensed at the BSM electrodes are often plotted 
as part of a map of the patient's chest. More particularly, 
the two-dimensional layout of the electrodes relative to 15 
the patient is first plotted. Then, for each point in the car- 
diac cycle at which measurements are made, the elec- 
trode sites having the same voltages are connected by 
contour lines, commonly referred to as "isopotential" 
lines. As a result, a series of isopotential maps are gen- 20 
erated for for each cardiac cycle. 

Although the generation of isopotential maps is rel- 
atively easy to accomplish, their interpretation can be 
difficult. In that regard, some relationship between the 
contours of the isopotential map and the particular car- 25 
diac features of interest must be established. Although a 
physician can make an initial evaluation based upon vis- 
ual comparisons with referential maps generated for 
populations of known physiology, the volume of data 
contained in the series of maps generated over one car- 30 
diac cycle can be overwhelming. Further, because map- 
ping is not presently widely performed, the physician's 
experience and familiarity with referential maps may be 
somewhat limited. 

A variety of techniques have been developed to 35 
reduce the complexity of map analysis. For example, 
one way of reducing the volume of data to be dealt with 
involves integration of the voltages sensed at each of 
the electrodes over an interval of time. In that regard, if 
the voltage at a given electrode is plotted as a function 40 
of time, the resultant waveform generally includes the P, 
Q, R, S, and T segments described above. By integrat- 
ing that voltage over, for example, the QRS or ST inter- 
val, a single value indicative of the sensed voltage at the 
electrode is produced. 45 

A single isointegral map, representative of cardiac 
activity over an entire cycle, can then be generated as 
follows. As with the isopotential map, the two-dimen- 
sional arrangement of the electrodes relative to the 
patient's chest is first plotted. Next, the integrals of the so 
waveforms at the various electrode sites are evaluated 
and those electrode sites having the same integral val- 
ues associated therewith are connected by lines, com- 
monly referred to as "isointegral" contours. As a result, 
a single isointegral map is produced representing car- 55 
diac activity over one cycle. 

Even with the data reduced to a single map per car- 
diac cycle, interpretation can be somewhat difficult. 
Each map includes a multitude of potentially significant 



contours to be evaluated by the physician. Trie experi- 
ence required to meaningfully interpret these contours 
based solely upon a visual review is considerably more 
extensive than that required to evaluate 12-lead ECGs. 
Such evaluation is further hampered by the relatively 
limited information available correlating map contours 
with cardiac activity. 

To enhance the value of such maps, increasing 
efforts are being made to identify significant map fea- 
tures and to develop systems that automatically extract 
such features from map data. With the appropriate fea- 
tures extracted, a physician can then more quickly and 
easily interpret the data. 

In that regard, one approach to feature extraction 
that is of particular interest is described in an article by 
Lux et al., entitled "Redundancy Reduction for Improved 
Display and Analysis of Body Surface Potential Maps, I. 
Spatial Compression," appearing in Volume 49 of Circu- 
lation Research at pages 186-196. Lux et al. disclose 
an approach in which a Karhunen-Loeve transform is 
used to compress the data to produce a basis function 
vector that allows the data to be reconstructed with min- 
imum representational error. A stepwise discriminant 
analysis is then applied to the vector coefficients to clas- 
sify the patient as diseased or normal. While the Lux et 
al. approach is a useful tool in the evaluation of BSM 
data, as will be described in greater detail below, the 
disclosed technique fails to consider the potential influ- 
ence of noncardiac patient variables on the analysis 
and does not adequately address the extraction of par- 
ticular features of interest from the data. 

In addition to reducing the complexity of BSM data 
analysis, attempts have been made to reduce the com- 
plexity of the data collection process. As will be appreci- 
ated, with 192 electrodes required to successfully 
generate a surface distribution of electrical activity, the 
time required to prepare a patient for data collection is 
significant. The electrode set required is further 
unwieldy, expensive, and may contribute to patient 
apprehension. Each of these factors tends to limit the 
utility of BSM data analysis. 

One work that is of particular interest in the area of 
"limited-lead" mapping is the article by Lux et al., enti- 
tled "Clinically Practical Lead Systems for Improved 
Electrocardiography: Comparison with Precordial Grids 
and Conventional Lead Systems," appearing in Volume 
59 of Circulation at pages 356-363. In this article, Lux et 
al. investigate the utility of several different lead sets in 
replicating the data collected with "conventional" 192 
electrode BSM systems. More particularly, sets of 32, 
30, and 9 leads are evaluated. Data collected from the 
limited leads is then transformed to represent data col- 
lected from a conventional lead set. Based upon the 
evaluation of this transformed data, Lux et al. conclude 
that 20-35 electrodes are required for BSM analyses. 

More particularly, Lux et al. emphasize numerous 
shortcomings of the nine lead data transformation. For 
example, Lux et al. note that the pattern error and root- 
mean-square (rms) error for the nine-lead set is larger 
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than the errors occurring when 32-lead arrays are used. 
Lux et al. conclude that the nine-lead array misses elec- 
trocardiographic information and that 20-35 electrodes 
are required for mapping. 

Efforts have also been devoted to the use of BSM 
data to identify particular cardiac anomalies of interest. 
For example, the subject of coronary artery narrowing is 
addressed in an article by Farr et al., entitled "Localiza- 
tion of Significant Coronary Arterial Narrowings Using 
Body Surface Potential Mapping During Exercise Stress 
Testing," appearing in Volume 59 of the American Jour- 
nal of Cardiology at pages 528-530, and an article by 
Sridharan et al., entitled "Use of Body Surface Maps to 
Identify Vessel Site of Coronary Occlusion," appearing 
in Volume 22S of the Journal of Electrocardiology at 
pages 72-81. As will be described in greater detail 
below, however, these approaches are limited by the 
form of the data analyzed, the feature extraction tech- 
niques used, and their ability to quantify, as well as 
localize, arterial narrowings. 

Another specific anomaly of interest is the heart's 
"inducibilrty", or susceptibility to externally induced 
arrhythmia. The subject of inducibilty is addressed to 
some extent , for example, in an article by Vincent et al. 
entitled "Use of QRST Area Distribution to Predict vul- 
nerability to Cardiac Death Following Myocardial Infarc- 
tion", appearing in Volume 68 of Circulation at page 
352, an article by Kubota et al. entitled "Relation of Car- 
diac Surface QRST Distribution to Ventricular Fibrilla- 
tion Thresholds in Dogs", appearing in Volume 78 of 
Circulation at pages 171-177, and an article by Han et 
al., entitled "Nonuniform Recovery of Excitability in Ven- 
tricular Muscle", appearing in Volume 1 4 of Circulation 
at pages 44-60. Again, however, the disclosed 
approaches are somewhat limited in nature and appli- 
cability. 

As will be appreciated from the preceding discus- 
sion, it would be desirable to provide a system that 
would allow a spatial distribution representative of, for 
example, the electrical activity of a heart to be analyzed 
without requiring the use of complex electrode configu- 
rations for data collection. The system should, however, 
take advantage of the added information content 
present in the spatial distribution. In addition, the sys- 
tem should format the spatially distributed data in a 
manner suitable for the particular feature to be 
extracted. The system should further be able to extract 
the features in a manner that distinguishes noncardiac 
influences on the data and that offers increased analyt- 
ical capabilities with respect to particular features of 
interest. 

Summary of the Invention 

In accordance with the invention, a system is dis- 
closed for determining whether a patient has a coronary 
disease, comprising: a first set of electrodes for receiv- 
ing a first set of cardiac signals from a patient, at less 
than twenty electrodes; a data conversion interface con- 



nected to receive the first set of cardiac signals and to 
convert the first set of cardiac signals into a digital form; 
memory means for storing a function that computes a 
canonical variable using three or more statistically 
5 determined coefficients taken from a set of coefficients; 
a central processing unit that is programmed to: 

a) transform the digitized set of cardiac signals into 
a greater number of digitized values, said set of val- 

10 ues representing a set of body surface mapping 
signals such that the first set of cardiac signals 
appear as if they had been received from the 
patient using a set of body surface mapping elec- 
trodes, at more than twenty electrodes, and said 

is body surface mapping electrodes had received the 
whole set of body surface mapping signals; said 
transformation using a mathematical expectation 
operator and stored referential data; 

20 b) to expand the transformed set of values repre- 
senting the of body surface mapping signals as a 
linear combination of predetermined spatial basis 
functions to determine said set of coefficients as 
being the coefficients of the linear expansion; 

25 

c) to compute a value for the canonical variable 
using the stored function; and 

d) to compare the value of the canonical variable to 
30 a cut point above which the canonical variable indi- 
cates that the patient has the coronary disease and 
below which the canonical variable indicates the 
patient does not have the disease; and 

35 a display that indicates whether the patient has the cor- 
onary disease, the display being activated by the central 
processing unit if the value of the value of the canonical 
variable is greater than the cut point. 

A method not in accordance with the invention is 

40 disclosed for monitoring and analyzing a first set of elec- 
trical signals, received from a patient at less than 20 
electrodes and conventionally analyzed on the basis of 
the signals' temporal characteristics. The method 
includes the step of transforming the first set of electri- 

45 cal signals into a transformed data set including infor- 
mation regarding the spatial distribution of the electrical 
signals relative to the patient. The method also includes 
the step of processing the transformed data set to eval- 
uate at least one patient feature of interest. 

so A method not in accordance with this invention is 
disclosed for processing information that is spatially 
representative of the electrical activity of a patient's 
heart to detect the presence of at least one cardiac 
characteristic of interest. The method includes the step 

55 of producing a data set that is usable in representing the 
information representative of the electrical activity of the 
patient's heart. A data subset is then identified from the 
data set for use in detecting the cardiac characteristics 
of interest. Finally, a feature value is computed from the 
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subset for use in detecting the presence of the charac- 
teristics. 

A method not in accordance with the invention is 
disclosed for reducing the influence of noncardiac 
parameters on data representing the spatial distribution s 
of electrical activity across a patient's chest. Another 
method is disclosed for evaluating the magnitude of a 
cardiac characteristic of interest. Additional methods 
are disclosed for detecting the presence of coronary 
artery stenosis in a patient, as well as the patient's sus- io 
ceptibilrty to externally induced arrhythmias. 

The invention is defined in the appended claims. 

Brief Descrip tion of the Drawings 

75 

The invention will presently be described in greater 
detail, by way of example, with reference to the accom- 
panying drawings, wherein: 

FIGURE 1 illustrates a 12-lead ECG/BSM system, 20 
constructed in accordance with this invention, 
attached to a patient; 

FIGURE 2 is a block diagram of the ECG/BSM sys- 
tem of FIGURE 1; 

FIGURE 3 is a flow chart illustrating the manner in 25 
which the system of FIGURE 1 receives data from 
the patient and transforms it into a conventional 
BSM format; 

FIGURE 4 is a graph of one lead of data, plotted as 
a function of time over one cardiac cycle; 30 
FIGURE 5 is a graph of the location of ten elec- 
trodes employed in the ECG/BSM system of FIG- 
URE 2, relative to a conventional 192-electrode 
mapping system; 

FIGURE 6 illustrates an isointegral map generated 35 
by the ECG/BSM system; 

FIGURE 7 is a flow chart illustrating one way in 
which the transformed BSM data produced in FIG- 
URE 3 is compressed for further analysis by the 
ECG/BSM system; 40 
FIGURE 8 is a flow chart illustrating one way in 
which clinical data is obtained from a patient popu- 
lation, including "normal" individuals, as well "dis- 
eased" individuals exhibiting the particular cardiac 
anomaly of interest; 45 
FIGURE 9 is a graph depicting the distribution of 
the patient population relative to the sensitivity and 
specificity of the analysis performed; 
FIGURE 1 0 is a graph of a receiver operating curve, 
providing an alternative depiction of the specifi- so 
city/sensitivity trade-off of FIGURE 9; and 
FIGURE 11 is a flow chart illustrating the manner in 
which the system of FIGURE 1 uses the information 
obtained in the analysis of FIGURE 8 to evaluate 
the data transformed and compressed in accord- 55 
ance with FIGURES 3 and 7. 



Detailed Description of the Preferred Embodiment 

A 12-lead ECG/BSM system 10, constructed in 
accordance with this invention, is shown attached to a 
patient 12 in FIGURE 1. The system 10 can be used to 
perform traditional 12-lead ECG analyses of the 
patient's heart. In addition, system 10 can be used to 
collect and analyze information regarding the spatial 
distribution of the heart's electrical activity. 

As will be described in greater detail below, the spa- 
tial analysis or body surface mapping performed by sys- 
tem 10 is accomplished with significantly fewer 
electrodes than was previously deemed possible. As a 
result, system 10 is simpler than conventional BSM sys- 
tems. Because system 10 can be attached to patient 12 
more quickly and easily than conventional BSM sys- 
tems, its clinical value is also enhanced. 

Trie manner in which the spatial data is processed 
by system 10 similarly offers several advantages over 
the feature extraction techniques conventionally 
employed. For example, by employing QRST integral 
information, the performance of a preferred arrange- 
ment of system 10 is not dependent upon the activation 
sequence exhibited by the patient's heart. Further, the 
system 10 may normalize the mapped data prior to fea- 
ture extraction to eliminate the effect of a variety of non- 
cardiac patient variables on the analysis. An additional 
enhancement involves the system's ability to detect, dis- 
tinguish, and quantify local cardiac anomalies. Yet 
another feature of the system 10 involves a technique 
that may be employed to reduce the complexity of the 
representations of the 'mapped data required for further 
analysis. 

Reviewing the construction of system 10 in greater 
detail, reference is had to the block diagram of FIGURE 
2. As shown, the ECG/BSM system 10 broadly includes 
an electrode cable set 1 4 and an ECG/BSM monitor 1 6. 
The electrode cable set 14 is attached to patient 12 and 
electrically connects the patient 12 to monitor 16. The 
monitor 16, in turn, senses the electrical activity of the 
patient's heart through the electrode cable set 14 and 
produces outputs that are used by a physician in making 
ECG-based and BSM-based diagnoses. 

The electrode cable set 1 4 is conventional in design 
and includes a four-electrode limb set 18 and a six-elec- 
trode precordial set 20. Addressing these components 
individually, the limb set 18 includes limb electrodes U 
R, F, and G, which are electrically connected to a limb 
set connector 22 by a four-conductor cable set 24. Elec- 
trodes L, R, F, and G are preferably designed to provide 
a low-impedance, adhesively maintained interface with 
patient 12. The connector 22 and cable set 24 provide 
electrically independent, conductive paths for the car- 
diac signals sensed at the four electrodes L, R, F, and 
G. 

The precordial set 20 includes six precordial elec- 
trodes p1, p2, p3, p4, p5, and p6. These six electrodes 
are connected to a common precordial set connector 26 
by a six-conductor cable set 28. The electrodes pi, p2, 
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p3, p4, p5, and p6 are also designed to provide a low- 
impedance, adhesively maintained interface with 
patient 12. The connector 26 and cable set 28 provide 
electrically independent, conductive paths for the car- 
diac signals sensed at the six precordial electrodes. As 5 
will be appreciated, it may be desirable to modify the 
number and location of electrodes used. 

Addressing now the construction of the monitor 1 6 
in greater detail, as shown in FIGURE 2, monitor 16 
includes a data conversion interface 30, central 10 
processing unit (CPU) 32 (such as an Intel 386 CPU), 
data storage system 34, input system 36, output system 
38, input interface 40 and output interface 42. Prefera- 
bly, these components are housed in a single unit. As 
will be appreciated, however, the monitor 16 may be a is 
system of discrete components including a personal 
computer, with supporting peripherals and interfaces. 

Addressing these various components individually, 
the data conversion interface 30 preferably includes a 
number of elements. For example, interface 30 may 20 
include an amplifier circuit 44, which amplifies the rela- 
tively low-level cardiac signals received from cable set 
14 to a level that is sufficient for subsequent processing. 
If desired, a filter circuit 46 having a bandwidth of 0.1 to 
250 Hertz, may also be included to remove noise from 25 
the received signals. Finally, an analog-to-digital (A/D) 
converter 48 is included to convert the analog cardiac 
signals into a digital form of suitable resolution for anal- 
ysis by the CPU 32. The interface 30 preferably includes 
nine channels to allow nine leads of data to be received so 
and processed simultaneously. 

The central processing unit (CPU) 32 is included to 
process the cardiac signals received from interface 30. 
As will be described in greater detail below, the CPU 32 
processes these signals with the aid of stored data and 35 
instructions, as well as various operator-generated 
inputs. The CPU 32 processes this information to gen- 
erate displays used by the physician in 12-lead ECG 
and BSM diagnosis. The data storage system 34 has a 
variety of responsibilities. In that regard, system 34 40 
retains the various program instructions that govern the 
manner in which CPU 32 processes the received car- 
diac signals. In addition, system 34 stores information 
about the cardiac signals in various formats as the sig- 
nals are being processed. Finally, the storage system 45 
34 stores referential data representing the cardiac activ- 
ity of normal and diseased patients for use by CPU 32 in 
interpreting the cardiac signals received from the elec- 
trode cable set 14. 

As will be appreciated, the data storage system 34 so 
may include random-access memory (RAM) and read- 
only memory (ROM) in chip form. The storage capability 
may be supplemented by the addition of floppy or hard- 
disk storage systems. Such components of the data 
storage system 34 may be included within, or external ss 
to, the monitor 16. 

The input system 36 is included to allow the opera- 
tor to provide data and instructions to the monitor 1 6. In 
that regard, the input system 36 may include a control 



panel or keyboard for use in controlling the monitor 16. 
Such a control panel allows, for example, the monitor 1 6 
. to be turned ON and OFF, ECG or BSM modes of oper- 
ation to be selected, a particular diagnostic approach to 
be chosen, or a particular display technique to be used. 
In addition, the input system 36 may be designed to 
allow program instructions to be written into storage 
system 34 through a keyboard or directly down-loaded 
from a remote storage system. 

The output system 38 includes a variety of devices 
useful to a physician in evaluating the performance of 
monitor 16 and the condition of the patient's heart. For 
example, the output system 38 may include displays 
that indicate whether the monitor 16 is ON and whether 
the ECG or BSM mode of operation has been selected. 
The output system 38 may further include a CRT or 
printer for displaying selected QRST waveforms used in 
ECG analyses, as well as isopotential and isointegral 
surface distributions used in BSM analyses. Finally, the 
output system 38 may include individual visible and 
audible displays indicative of particular cardiac condi- 
tions detected by the CPU 32, as described in greater 
detail below. 

The input interface 40 connects the input system 36 
to the CPU 32. As will be appreciated, interface 40 is 
designed to condition input signals so that they can be 
applied to, and properly interpreted by, CPU 32. The 
output interface 42, in turn, connects the output system 
38 to CPU 32 and is designed to condition the outputs 
of CPU 32 so that they are usable by the output system 
38. 

Having reviewed the basic components of the 
ECG/BSM system 10. the operation of system 10 will 
now be considered in greater detail. In that regard, the 
electrode cable set 14 is attached to patient 12 in 
accordance with conventional 12-lead protocols. As 
previously described in detail, the limb electrodes L, R, 
F, and G are connected to the patient's limbs and the 
precordial electrodes p1 , p2, p3, p4, p5, and p6 are con- 
nected to the patient's chest. The positioning of the 
electrodes influences the strength of the received sig- 
nals and, thus, is important in ensuring the proper anal- 
ysis of cardiac activity. 

With the electrode cable set 14 properly attached to 
patient 12, the monitor 1 6 offers several different modes 
in which the patient's cardiac activity can be evaluated. 
For example, a physician can instruct the monitor 16 to 
operate in an ECG mode by appropriately actuating 
controls in the input system 36. In this ECG mode, the 
monitor 16 causes the output system 38 to selectively 
display the twelve leads of ECG data for analysis by the 
physician. The operation of system 10 in the ECG mode 
is conventional and need not be discussed in further 
detail here. 

The operation of system 10 in a BSM mode is, how- 
ever, of particular interest and is described in greater 
detail in connection with the flow chart shown in FIG- 
URE 3. In that regard, the physician first instructs the 
monitor 16 to operate in the BSM mode by actuating the 
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appropriate control in the input system 36, as shown in 
block 50 of the flow chart. In the BSM mode, monitor 16 
receives analog signals representative of the potentials 
at the L, R, F, p1 , p2, p3, p4, p5, and p6 electrodes, 
measured relative to the potential of the ground elec- 5 
trode G. 

These nine signals are initially applied to the data 
conversion interface 30 at block 52. As noted previously, 
the interface 30 is preferably a nine-channel system, 
allowing each of the signals to be continuously received 10 
and processed by the monitor 16. The amplifier circuit 
44 amplifies the signals at block 54 and the filter circuit 
46 removes undesirable frequency components there- 
from at block 56. Finally, at block 58. the A/D converter 
48 converts the analog signals into digital form. 15 

From this point on, the operation of system 10 
depicted in the flowchart of FIGURE 3 represents steps 
performed by CPU 32 in response to program instruc- 
tions stored in system 34. The first of these operations, 
shown in block 60 is the removal of offset components 20 
to establish a baseline from which signals can be more 
readily compared and analyzed. Next, at block 62, each 
of the nine signals are simultaneously sampled and 
stored in system 34. at the conventional 250 Hertz sam- 
ple rate used by 12-lead ECG systems. 25 

The next series of preprocessing steps performed 
by the CPU 32 is designed to identify when the received 
cardiac signals are suitable for further analysis. At block 
64. the CPU 32 evaluates, for example, one of the nine 
signals over a single cardiac cycle. The CPU 32 gener- 30 
ates a waveform template having certain amplitude and 
temporal characteristics that are indicative of the signal 
selected, during the cycle reviewed. 

Then, at block 66, the CPU 32 compares the signal 
received over several consecutive cardiac cycles with 35 
the template, if the signal continues to exhibit a periodic 
waveform whose amplitude and temporal characteris- 
tics are within set tolerances of the template values, the 
signal is stable and is considered to be an accurate rep- 
resentation of cardiac activity. As a result, the CPU 32 40 
will be directed to analyze all nine signals collected over 
the last cardiac cycle, at block 68. 

On the other hand, if the template parameters are 
not satisfied by the signal over the consecutive cycles 
reviewed at block 66, the CPU 32 is directed to restore 45 
operation to block 64. As a result, a new set of template 
parameters will be generated at block 64. The signal will 
then be compared to the new template over several sub- 
sequent cardiac cycles at block 66. Eventually, the tem- 
plate parameters will be satisfied and the operation of so 
CPU 32 will proceed to block 70. 

At block 70, CPU 32 reviews the nine signals col- 
lected during the cardiac cycle selected by the template 
comparison process reviewed above. FIGURE 4 is an 
illustration of one of the nine signals, for example, the 55 
potential at electrode p1 relative to electrode G, plotted 
over one cardiac cycle. As shown, the signal exhibits 
the conventional PQRST segments. 

In further processing the nine signals, the CPU 32 
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first computes the twelve leads of data used in 12-lead 
ECG analysis. Thus, the conventional leads are readily 
available for display when the ECG/BSM system 10 is 
used in the ECG mode. 

The CPU 32 also identifies the portion of the nine 
signals to be analyzed by system 1 0 in the BSM mode. 
Assuming first that the QRST portion of each signal is to 
be evaluated, the CPU 32 identifies the beginning of the 
Q segment and the end of the T segment of each signal. 
The identification of these points can be accomplished 
by monitoring the signals for slope reversals occurring 
at intervals in which such transitions are expected to 
occur. With the beginning of the Q segment and the end 
of the T segment established as end points, the CPU 32 
then limits further analysis of the nine signals to those 
portions of the signals occurring between the end 
points. 

As will be appreciated, the QRST intervals of the 
conventional twelve leads could be computed in place 
of the corresponding intervals of the nine electrode sig- 
nals. Because the desired spatial information is con- 
tained directly in the nine signals, however, the 
complexity of subsequent processing is reduced by 
eliminating the twelve-lead conversion. 

Once the QRST intervals of interest have been 
determined, the CPU 32 proceeds to integrate each of 
the nine intervals at block 72. As will be appreciated, 
nine numeric values, corresponding to the area under 
the nine QRST curves, are thus computed for the car- 
diac cycle being reviewed. As a result, the complexity of 
the information to be further processed is reduced by a 
factor of 250. 

Another advantage of using the QRST area, as 
opposed to the waveform itself or some limited portion 
of the waveform, is that the QRST area is relatively 
independent of the activation sequence of the myocar- 
dial cells. Thus, an anomaly in myocardial conduction, 
such as, for example, a bundle brand blockage, that is 
unrelated to the cardiac feature of interest will not alter 
the QRST area computation or any further analysis 
based upon that computation. 

At block 74, the CPU 32 is called upon to begin a 
transformation of the nine values computed at block 72 
for use in preparing a body surface map. In that regard, 
the essence of this transformation is the conversion of 
data that is conventionally analyzed on a temporal basis 
into a format for use in a spatial analysis. 

By way of illustration, assume that a 1 92-electrode 
set represents the conventional standard number of 
electrodes required to fully evaluate the spatial distribu- 
tion of signals across the patient's chest. The pattern of 
these 192 electrodes is shown graphically in the map of 
FIGURE 5 and includes 16 columns of 12 sites each. 

Of these 192 "standard" electrode sites, the nine 
sites whose physical location corresponds most closely 
to the location of properly placed conventional elec- 
trodes L, R, F, p1, p2, p3, p4, p5, are identified. In FIG- 
URE 5, these sites are identified with the letter x. With 
the appropriate limited-lead set identified, a linear trans- 
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formation of data from the nine electrode sites to 192 
electrode sites can be accomplished in the manner 
taught by Lux et al., cited above. 

More particularly, at block 74 the CPU 32 computes 
a transformed vector A1, which includes the set of 192 
transformed values, as follows: 

A1=E(A2)+(T1)(A3-(E(A3))/T2 ; where E is a mathe- 
matical expectation operator; A2 is a vector including a 
referential 192 lead data set stored for use in performing 
the transformation; A3 is a vector including the nine val- 
ues measured for the patient; T1 is the cross-covari- 
ance matrix of A2 and A3; and T2 is the covariance 
matrix of A3. 

On a functional level, this process involves a linear 
transformation employing a 9-row, 183-column matrix 
retained in storage system 34. This transformation 
matrix is computed from data measured with a 192- 
electrode set and each point in the matrix is a least- 
squares approximation of the value required to ensure 
that the sum of the possible squared prediction errors, 
occurring when the transformation is performed, is kept 
to a minimum. 

Once the transformed vector has been determined, 
CPU 32 preferably normalizes the vector to unit length 
at block 76. As will be appreciated, the magnitudes of 
the electrical signals received at the various electrodes 
are a function of a variety of factors. Some of these fac- 
tors are related to the operation of the patient's heart 
and may be indicative of a particular anomaly of inter- 
est. On the other hand, other factors, including such vol- 
ume-conductor related parameters as the patient's 
weight, age, and sex, may also influence the amplitude 
of the signals received. As a result, the signals received 
from two patient's with identical heart conditions may 
exhibit significantly different amplitudes. 

The normalization step removes this form of non- 
cardiac differentiation between patients. In that regard, 
each element of the 192 element vector is squared and 
their sums added, with the magnitude of the vector 
being equal to the square root of that quantity. By divid- 
ing each of the 192 elements by the vector magnitude, 
the vector is then normalized to have a magnitude of 
one. As a result, the patterns of the signals are retained 
but the potentially skewing influence of noncardiac 
patient parameters is eliminated. As will be appreciated, 
information regarding the magnitude of the signal can 
still be retained, if desired. 

Once the normalized, transformed vector has been 
determined, CPU 32 may instruct the output system 38 
to display the vector as a series of isointegral lines at 
block 78. In that regard, the transformed vector includes 
192 numbers that are representative of the QRST area 
that would have been expected to be measured if 192 
electrodes had been employed in place of the nine 
relied upon. The output system 38 displays a map, 
shown in FIGURE 6, representative of the 192 electrode 
sites and joins those sites having substantially the same 
number assigned thereto by lines. As a result, the map 
of the patient's chest includes a set of isointegral con- 
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tour lines indicative of the integrals of the QRST areas 
that would be expected to be measured at points along 
the contour lines. 

A physician who is experienced in reviewing the 

5 body surface maps of patients having known cardiac 
diseases can visually analyze the mapped transforma- 
tion data to look for features potentially indicative of car- 
diac disease. As noted previously, however, the work by 
Lux et al. indicates that such analyses would not be 

w useful because the correlation between data trans- 
formed from nine leads and data measured with 192 
leads is not sufficiently high. However, despite the low 
correlation, applicants have discovered that nine-lead 
transformed data can be quite useful in detecting and 

15 analyzing certain cardiac anomalies. 

The preceding discussion describes one approach 
to the collection of information concerning the spatial 
distribution of electrical activity across the patient's 
chest. By successfully using an electrode configuration 

20 that is conventionally used to provide only a temporal 
evaluation of the electrical activity, the physical and 
operational complexity of the spatial analysis is advan- 
tageously reduced. The following discussion of the 
extraction of features from the transformed BSM data 

25 applies equally well, however, to the analysis of data 
collected from conventional BSM electrode sets. 

Addressing now the manner in which the analysis 
of the transformed BSM data by ECG/BSM system 10 
proceeds, reference is had to FIGURE 7. Although the 

30 physician can simply analyze the data based upon a 
visual review, the extent and subtleties of the informa- 
tion content in the map, along with the limited experi- 
ence of most physicians in the review of BSM data, 
make it desirable for the CPU 32 to perform additional 

35 feature extraction and analyses. 

To that end, the CPU 32 first attempts to further 
reduce the volume of information to a more manageable 
level. One such approach involves the representation of 
the mapped data by a set of basis functions, as indi- 

40 cated at block 80. These basis functions are essentially 
a set of building blocks that can be used to reconstruct 
the mapped data. The advantage of using the basis 
functions is that the functions, as well as the way in 
which they must be combined to reconstruct the data 

45 set, can be stored and evaluated more easily than the 
mapped data itself. 

The most commonly used basis functions in non- 
mapping data compression applications are sinusoidal 
functions used to reconstruct data on the basis of its fre- 

so quency content. However, in many instances, it is possi- 
ble to identify nonsinusoidal basis functions that are 
more useful either because they are computationally 
more efficient or because they permit greater data com- 
pression. 

55 In one embodiment, the CPU 32 identifies the most 
representative basis functions for subsequent analysis 
at block 80. In that regard, a Karhunen-Loeve (KL) 
expansion is used to identify those basis functions that 
allow the difference between the mapped data and a 
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reconstruction achieved with the basis functions, i.e., 
the representation error, to be minimized. The basis 
functions selected are the eigenvectors of the covari- 
ance matrix of the data, ordered by decreasing eigen- 
value magnitude. KL expansions are described in 5 
greater detail in, for example, Chen, Statistical Pattern 
Recognition, pages 36-38 and 161. 

The KL expansion results in the production of 192 
basis functions. In the preferred arrangement, however, 
only the twelve highest order functions are used at block 10 
82. The original 192-dimension map M can then be 
expressed as: 

12 

Z ci 

1=1 15 

where <j>i represents the basis functions selected and ci 
represents coefficients that can be used with the basis 
functions to reconstruct the map M. 

These twelve basis functions are small enough in 
number to be manageable from a processing stand- 20 
point, while still preserving a sufficient amount of infor- 
mation to be useful in evaluating the particular cardiac 
anomaly or anomalies of interest. Although error theo- 
retically decreases with the inclusion of additional basis 
functions, experience indicates that the response 25 
becomes relatively flat when more than ten functions 
are used and, for example, the use of 1 8 basis functions 
does not appear to provide any better results than 12 
basis functions. 

Although the use of the KL expansion to minimize 30 
representation error, as described above, is valuable, it 
may not be the optimal approach. In that regard, the fact 
that this expansion allows the most accurate represen- 
tation or reconstruction of the mapped data to be 
achieved does not mean that it is also always the most 35 
useful expansion for classifying cardiac anomalies. For 
example, a different expansion could be selected that 
might increase the representation error but would allow 
a particular cardiac anomaly of interest to be more read- 
ily classified. 40 

In that regard, one modification of the expansion 
that may be particularly suitable for use in this applica- 
tion is that described by Fukunaga et al. in an article 
entitled "Application of K-L Expansion to Feature Selec- 
tion and Ordering: A Criterion and an Algorithm for 45 
Grouping Data," appearing in Volume C19 of IEEE 
Transactions on Computers at page 311. The disclosed 
approach is a modification of the conventional KL 
expansion involving the transformation of the covari- 
ance matrix before the eigenvectors are extracted. so 

The basis functions identified by this modified KL 
expansion are selected to enhance the differences 
between two pattern classes. More particularly, in the 
present application, they are selected to accentuate the 
difference between a "normal" class that does not 55 
exhibit the cardiac anomaly of interest and a "diseased" 
class that does exhibit the anomaly. The basis functions 
identified may not, however, be the best at reconstruct- 
ing the mapped data itself. 



As noted previously, however, in the current 
arrangement the conventional KL expansion is used to 
achieve the best representation of the mapped data. 
This approach allows the ECG/BSM system 10 to 
remain relatively adaptable because the initial emphasis 
on representation is believed to retain information that 
might be useful in discriminating various features. On 
the other hand, if the expansion were selected to 
enhance the classification of a particular anomaly, the 
utility of the expansion in subsequently classifying other 
anomalies could be compromised. 

With the appropriate twelve basis functions estab- 
lished at block 82, the transformed integral data is then 
multiplied by the orthonormal basis functions at block 84 
to produce twelve coefficients. The analysis next typi- 
cally involves some form of statistical treatment of the 
coefficients. More particularly, even though significant 
processing has been performed to reduce the data pre- 
sented for analysis, it may still be of sufficient scope and 
subtlety to render meaningful analysis difficult. As a 
result, some form of stepwise discriminant analysis, 
logistic regression, or cluster analysis, discussed gener- 
ally in AfHi et al., Computer-Aided Multivariate Analysis, 
parts II and III, may be performed. 

In that regard, the subsequent statistical analysis of 
the twelve coefficients is based upon information 
obtained from "normal" and "diseased" patient popula- 
tions in a clinical setting. This information is stored in the 
data storage system 34 for use by the ECG/BSM sys- 
tem 10. 

Reviewing the manner in which this information is 
obtained in greater detail, reference is had to FIGURE 8 
in which the various steps to be performed are illus- 
trated in flow chart form. First, suitable normal and dis- 
eased populations are selected as indicated at block 86. 

The normal population includes a statistically ade- 
quate number of patients that are in good health, have 
had no history of heart disease, and have been evalu- 
ated for the absence of a particular cardiac anomaly or 
anomalies of interest by a generally accepted clinical 
technique or techniques. The diseased population, on 
the other hand, includes a statistically adequate number 
of patients that have been identified as exhibiting the 
particular anomaly or anomalies of interest. Where 
appropriate, a quantification of the clinically evaluated 
anomaly exhibited by each patient in the diseased pop- 
ulation is also determined. 

With suitable populations identified, data represent- 
ing the electrical activity of each patient's heart is col- 
lected at block 88. In that regard, the same electrode 
configuration is employed and the signals sensed at the 
electrodes are collected in the same manner as those 
that will eventually be evaluated by the ECG/BSM sys- 
tem 10. Similarly, at block 90. the signals are processed 
in the manner described above in connection with 
blocks 52 through 84. 

Thus, the collection and processing of data from the 
normal and diseased populations represented by blocks 
88 and 90 can be performed by a prototype ECG/BSM 
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system 10. For the arrangement described above, the 
data collected from each patient in the normal and dis- 
eased populations will be reduced to a set of twelve 
coefficients. 

Next, the twelve coefficients associated with each 5 
of the patients in both populations, as well as any perti- 
nent clinical data concerning the patients, are down- 
loaded from the data storage system 34 of the prototype 
ECG/BSM system 10 to a memory system used by a 
suitable auxiliary computer at block 92. The computer is 10 
programmed with, for example, the BMDP LR (logistics) 
or TM (linear) software packages available from Statisti- 
cal Software of Los Angeles, California and performs a 
statistical analysis of this data, as described in greater 
detail below, to identify the particular coefficients that is 
are most useful in detecting the particular cardiac 
anomaly or anomalies of interest. 

More particularly, the mainframe computer per- 
forms a stepwise regression analysis on the coefficients 
stored for all of the patients in the normal and diseased 20 
populations as indicated at block 94. A separate such 
regression analysis is performed to for each feature of 
interest. As an initial step in the analysis, the computer 
first looks for the single one of the twelve coefficients 
that best allows the presence or absence (and perhaps 25 
extent) of the particular anomaly of interest to be identi- 
fied for all individuals in the population. 

In that regard, a canonical variable V or feature 
value is computed for each patient in the population. 
The canonical variable V is used to determine the ability 30 
of the coefficient(s) selected to accurately allow dis- 
eased and normal patients to be distinguished. If the 
computer employs a linear analysis, the variable V is 
calculated as: 

35 

w0(c0)+w1 (c1)+...+wn(cn) 

where cO through cn are the coefficients identified by 
the stepwise analysis and wO through wn are weighting 
factors associated with each coefficient. Note that all but 40 
one of the coefficients is set equal to zero during this ini- 
tial phase of the discriminant analysis. Also, if a logistic 
analysis is employed, the canonical variable V is used in 
the logarithmic analysis: 

e v /(l + e v ) 46 

As part of the analysis of the canonical variable, the 
computer also makes assumptions about a decision 
boundary that divides the normal and diseased patient so 
populations. If the canonical variable computed for a 
particular member of the population and particular coef- 
ficient is greater than the decision boundary or outpoint, 
then the patient is considered to exhibit the cardiac 
anomaly of interest. On the other hand, if the variable is 55 
less than or equal to the outpoint, then the patient is 
considered normal. The number of incorrect diagnoses 
made by this process during the evaluation of the popu- 
lation is then counted to estimate the probability of mis- 
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take for the particular combination of cutpoint and 
coefficient selected. 

Ultimately, the computer identifies the single coeffi- 
cient cO that will minimize the likelihood of incorrectly 
distinguishing diseased and normal patients in the clini- 
cal population for all possible outpoints. As part of the 
stepwise nature of the analysis, the computer next 
determines whether the predictive ability of the analysis 
can be enhanced by including another one of the coeffi- 
cients in the computation of the canonical variable V. If 
so, the "best" additional coefficient is selected as a sec- 
ond coefficient to be used by system 10. This process is 
repeated until no additional coefficients can be used to 
enhance the detection and analysis of the anomaly of 
interest. Ultimately, the mainframe computer may deter- 
mine that, for example, the third, fifth, and seventh of the 
twelve coefficients are the best for use in evaluating the 
anomaly of interest. 

Once the appropriate coefficients are identified by 
the computer, their identification is programmed into the 
storage system 34 of the ECG/BSM system 10 for use 
in analyzing the transformed patient data and providing 
a suitable diagnostic output to the operator. Information 
concerning the nature and magnitude of the anomaly 
associated with the coefficients is also stored. If several 
different anomalies are to be evaluated, the coefficients 
to be used in computing each of the different canonical 
variables (V1 , V2, through Vn) are stored, along with the 
related information concerning the nature of the anom- 
alies. 

Before leaving this discussion of the interpretation 
of the patient population data, the results of the popula- 
tions analysis will be graphically recapped. In that 
regard, FIGURE 9 illustrates a plot of the distribution of 
the canonical variables V computed for all of the 
patients in the normal and diseased populations using 
the coefficients selected by the mainframe computer as 
part of the stepwise regression analysis. The extreme 
left-hand side of the curve indicates maximum sensitiv- 
ity (i.e., sensitive because all patients have a canonical 
variable greater than zero) and zero specificity (i.e., not 
specific because both normal and diseased patients are 
equally likely to have a variable V greater than zero). On 
the other hand, the extreme right-hand side of the curve 
indicates maximum specificity (i.e., only truly diseased 
patients will have a variable V equal to one) and zero 
sensitivity (i.e.. no patients will actually have a variable 
equal to one). 

As represented by the solid curve A in FIGURE 9. 
the distribution of variables for a linear analysis will typ- 
ically include two spaced-apart. bell-shaped sections, 
which generally separately represent the normal and 
diseased populations. By way of illustration, the varia- 
bles associated with the normal population may be cen- 
tered, for example, around 0.2, while the variables 
associated with the diseased population are centered 
about 0.8. A "outpoint" can be established between 
these two sections of curve A at, for example, 0.5, with 
variables greater than 0.5 being associated with a 
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patient exhibiting the anomaly of interest and variables 
less than or equal to 0.5 being associated with a normal 
patient. If a logistic analysis is used in place of the linear 
analysis, a favorable spreading of the normal and dis- 
eased sections of the curve with respect to the outpoint 5 
may be achieved, as indicated by curve B of FIGURE 9. 

As will be described in greater detail below, the cut- 
point may be of use to the ECG/BSM system 10 in its 
analysis. In that regard, the system 1 0 may compute the 
canonical variable for a given patient, using the particu- 10 
lar ones of the patient's coefficients identified by the dis- 
criminant analysis. A threshold corresponding to the 
outpoint can then be used to interpret whether the 
canonical variable V indicates the presence of the par- 
ticular anomaly of interest. Patients having a variable V 75 
greater than the threshold can be diagnosed as suffer- 
ing from the anomaly of interest, while the remaining 
patients are diagnosed as normal. 

By adjusting the outpoint and threshold to, for 
example, 0.3 the sensitivity of the analysis can be 20 
increased. This modification might be appropriate when 
the ECG/BSM sytem 10 is used as an initial screening 
tool to identify potentially diseased patients for further 
testing because the operator would likely want to reduce 
the possibility of a diseased patient being diagnosed as 25 
normal. On the other hand, if the ECG/BSM system 10 
is being used to identify patients needing surgery, the 
outpoint and threshold used might be adjusted upward 
to, for example, 0.8 to increase the specificity of the 
analysis. Thus, the possibility of a normal patient being 30 
subjected to surgery would be reduced. 

An alternative way of depicting the tradeoff between 
sensitivity and specificity inherent in the analysis is illus- 
trated in FIGURE 10. In that regard, FIGURE 10 illus- 
trates a receiver operating curve (ROC). The ROC is a 35 
plot of the data from FIGURE 9, with the y-axis repre- 
senting sensitivity, or true positive fraction, and the x- 
axis representing the quantity of M 1 -specificity," or false 
positive fraction, both expressed as percentages. 

The area under the ROC curve of FIGURE 10 is 40 
representative of the performance of the ECG/BSM sys- 
tem 10 with respect to the identification of the particular 
anomaly of interest. In that regard, if the particular coef- 
ficients used by the ECG/BSM system 10 could allow 
the system 10 to correctly diagnose all patients, the 4s 
analysis would have a specificity of one and a sensitivity 
of one. The area under the curve of FIGURE 10 would 
then be equal to one. Thus, the area under the curve 
generated for a real system will be some fraction of one. 

Returning now to the detailed discussion of the so 
analysis of an individual patient by the ECG/BSM sys- 
tem 10, as will be recalled from block 84 of FIGURE 7, 
the system 10 computes twelve coefficients associated 
with the electrical activity of the patient's heart. An indi- 
cation of the particular ones of these coefficients (e.g., 55 
the third, fifth, and seventh) to be used in the analysis of 
an anomaly of interest has been stored in the data stor- 
age system 34 of system 1 0 at block 96 of FIGURE 8. In 
addition, information concerning various outpoints suita- 



ble for use by system 10 in different applications has 
been stored, along with information concerning the 
nature and extent of anomalies associated with the var- 
ious canonical variables to be computed. 

Referring now to the flow chart of FIGURE 1 1, the 
ECG/BSM system 10 uses the stored information to 
determine the canonical variable V1 associated with the 
anomaly of interest at block 100. In that regard, if the 
analysis performed by the auxiliary computer in FIG- 
URE 8 indicated that the third, fifth, and seventh coeffi- 
cients allow the anomaly to be detected with the lowest 
likelihood of error, then the ECG/BSM system 10 com- 
putes the variable V1 using the third, fifth, and seventh 
coefficients developed from the patient's data. Specifi- 
cally, those coefficients are substituted into the appro- 
priate linear or logarithmic formula for the canonical 
variable V. 

At block 1 02, the canonical variable V1 is then com- 
pared against an appropriate outpoint, which may be 
input by the operator either manually or as part of a 
mode selection process. If the outpoint is exceeded, the 
patient is deemed to exhibit the anomaly of interest. On 
the other hand if the variable is less than or equal to the 
outpoint, then the patient is considered normal. 

As will be appreciated, the statistical analysis per- 
formed by the auxiliary computer included a considera- 
tion of the outpoints to be used and their influence on 
the performance of the analysis. Thus, a table listing the 
success of the coefficients selected at correctly identify- 
ing the anomaly for the full range of outpoints selectable 
can be provided to data storage 34 for use in the analy- 
sis. In that regard, the CPU 32 might be instructed to 
limit the range of outpoints selectable to those that will 
ensure a certain predictive accuracy for the coefficients 
being used. 

As a more sophisticated enhancement, storage 34 
might be provided with multiple sets of coefficients that 
provide the best predictive ability for different cutpoints. 
The CPU 32 can then be programmed to review the cut- 
point selected and decide whether an alternative set of 
coefficients could be to used to achieve a greater accu- 
racy. 

As will be appreciated, the resultant NORMAL ver- 
sus DISEASED diagnosis can be displayed by the out- 
put system 38 of the ECG/BSM system 10 at block 104 
of FIGURE 11 . In addition, display of the numeric value 
of the canonical variable V1 and the outpoint can be 
provided to enhance the operator's understanding of the 
evaluation performed. Further insight into the analysis 
performed may also be achieved by plotting the 
patient's variable V1 on a display of the curve of FIG- 
URE 9 or 10. 

The ECG/BSM system 10 can also be programmed 
to combine different groups of the patient's coefficients 
to compute different canonical variables, such as V2 
and V3, associated with different anomalies of interest 
in the same patient. The output system 38 would then 
include suitable displays of each evaluation. 

Jn addition to the basic analysis described above, 
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clinical experience suggests that the assessment of the 
canonical variable V made by the ECG/BSM system 10 
can be further enhanced by the inclusion of certain clin- 
ical factors that are not present in the electrical signals 
sensed at the electrode sites. For example, if the analy- 5 
sis of a patient's cardiac activity in the foregoing manner 
indicates a low probability of ischemia but the patient is 
experiencing chest pain, there is perhaps a greater like- 
lihood of the probability being inaccurate than if no 
chest pain is experienced. The same is true for a patient w 
having a personal history of myocardial infarction or a 
family history of cardiac disease. 

The ECG/BSM system 10 can be programmed to 
account for these clinical factors, as well as others such 
as sex, age, and cholesterol levels. Specifically, an addi- is 
tional term corresponding to each clinical factor to be 
considered may be added to the computation of the 
canonical variable V. Alternatively, the threshold level 
against which the variable V is compared for diagnostic 
purposes may be adjusted as a function of the clinical 20 
factors present. It is expected that the inclusion of these 
clinical factors may result in an adjustment of up to 30 
percent in the variable or cutpoint used. 

The basic analysis described above for extracting 
features from the mapped data can be altered in various 25 
ways. One important variation involves the elimination 
of the transformation step. In that regard, it is possible to 
obtain a suitable spatial distribution of charge relying 
simply upon the twelve lead data collected. As will be 
appreciated, this simply requires a modification of the 30 
post-transformation processing in accordance with the 
general principles outlined above. 

The inclusion of software to transform the twelve- 
lead data to 192-lead data does, however, offer an 
advantage. Specifically, alternative electrode sets 35 
including more or less than the ten electrodes dis- 
cussed above can be used and the transformation proc- 
ess will cause the data to be converted to a form that is 
compatible with the existing post-transformation 
processing software. On the other hand, if no transfer- 40 
mation were performed, alternative processing would 
be required to allow the ECG/BSM system 10 to be 
used with electrode sets having more or less than the 
ten electrodes. 

Another important variation in the analysis per- 45 
formed by the ECG/BSM system 10 relates to the 
amount, or form, of the electrical activity sensed at the 
electrodes. In that regard, although the analysis 
described above is based upon the evaluation of QRST 
integrals, as will be appreciated, various other formats so 
can be analyzed depending, in part, upon the nature of 
the cardiac anomalies of interest. For example, in addi- 
tion to the QRST interval, the QRS and ST intervals can 
be analyzed. Also, the data collected during these inter- 
vals can be analyzed on a point-by-point basis or as a ss 
single integral representative of the entire interval. 

As noted above, one determinant for the format of 
the signals to be analyzed is the type of cardiac charac- 
teristic to be evaluated. In that regard, as was previously 
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described, the features of the electrical waveform 
sensed at any given electrode are directly related to the 
physiological activity of the heart. Thus, it is reasonable 
to expect that the particular segment of the waveform to 
be analyzed, as well as the way in which the segment is 
represented, would involve the identification of the por- 
tion of the waveform that is physiologically most closely 
related to the particular anomaly of interest. 

As one illustration of the relationship between the 
waveform and the activity to be analyzed, many forms of 
the cardiac disease to be analyzed can be grouped as 
either activation-related diseases or recovery-related 
diseases. Activation-related diseases include, for exam- 
ple, necrosis, accessory pathway blockages, and bun- 
dle blockages. Recovery-related diseases, on the other 
hand, include ischemia and nerve disorders. Generally, 
activation-related diseases are considered to have the 
greatest impact upon the QRS segment of cardiac 
waveforms, while the ST segment is generally consid- 
ered to be most representative of recovery-related dis- 
eases. Because the processes of activation and 
recovery are related, the entire QRST interval is some- 
times also considered to be representative of recovery- 
related features, especially those involving the disparity 
of ventricular recovery properties. 

The analysis of the integral or area of the QRST 
interval is reviewed above. If the data is not compressed 
by integration, the analysis can still proceed in much 
much the same manner. For example, instead of pro- 
ducing a single isointegral map, the ECG/BSM system 
10 may generate a plurality of isopotential maps corre- 
sponding to each of the discrete segments of the inter- 
val. Basis functions may also be identified for for each 
discrete segment of the interval and analyzed through 
the use of a similar statistically constructed software 
package stored in the memory 34 of ECG/BSM system 
10. 

If, for example, each 250 to 400 segments of the 
QRST interval are analyzed individually, rather than as 
a single integral, the complexity of the processing per- 
formed by the ECG/BSM system 10, as well as the sta- 
tistical analysis required to program system 10, 
increases by an order of magnitude. However, the 
advantage of this approach over the use of integrals is 
that the analysis is now abte to account for temporally 
specific information in the waveform. 

The analysis of the QRS interval parallels the anal- 
ysis of the QRST interval, described above and may be 
accomplished in several different ways. For example, 
the map of the QRS integral, or the maps of the individ- 
ual segments of the QRS interval, may be displayed at 
the output system 3d and visually analyzed by an expe- 
rienced physician. As will be appreciated, the CPU 32 
can be easily programmed to operate in the same man- 
ner described above, except that the T segments of the 
waveforms are eliminated prior to further processing 
and analysis. This analysis may be useful in identifying, 
for example, the presence of necrotic tissue, ventricular 
hypertrophy, and previous myocardial infarction. 
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As represented by block 104 of FIGURE 11, ven- 
tricular activation time (VAT) maps can be generated 
and evaluated for the QRS interval. In that regard, the 
QRS portion of the waveform is most sensitive to ven- 
tricular activation. The VAT map provides a useful way 5 
of comparing the earliest and latest parts of the ventricle 
to activate, which may be indicative of the presence of, 
for example, ischemia or beat irregularities. 

A VAT map is generated by using the fully 
expanded waveform data, prior to the performance of a 10 
basis function compression at block 78. The maximum 
negative derivatives of the waveforms corresponding to 
the various electrode sites are then computed. As will 
be appreciated, these derivatives represent ventricular 
activation. By displaying the maximum derivatives as 75 
isoderivative contours, a map of the earliest and latest 
portions of ventricle to activate is provided. 

When the ST segments are to be evaluated, the 
CPU 32 must be programmed to alter the portions of the 
waveforms used in the analysis. The processing of the 20 
ST segment integral or waveform data then proceeds in 
the manner described above for the QRST and QRS 
intervals. One feature of particular interest that may be 
suitably identified by ST segment analysis is the detec- 
tion of cardiac injury currents. In that regard, it has been 25 
found that, for example, disease-related injury or trauma 
to cardiac tissue may introduce a displacement or shift 
in the ST segment. This shift is most readily observed 
by a comparison of the plots of isopotential contour 
maps, which will readily illustrate temporal variations in 30 
the activity of the heart. 

One other approach used to evaluate the ST seg- 
ments involves the use of snapshots of isopotential 
frames. As will be appreciated, this technique requires 
the CPU 32 to retain the waveform data from the various 35 
electrodes collected over the cardiac cycle, before com- 
pression by the computation of integrals. Isopotential 
maps are then generated at block 78 in FIGURE 3 for 
certain predetermined times during the ST segment of 
the cardiac cycle, such as 40, 60, or 80 milliseconds 40 
after the initiation of the interval. As will be appreciated, 
these snapshots can be used to evaluate injury-related 
ST segment displacements, as described above. 

Turning now to several specific applications of inter- 
est, one cardiac anomaly that has been found to be par- 45 
ticularly suitable for detection, using 12-lead 
transformed BSM data and the feature extraction tech- 
niques described above, is coronary artery stenosis. As 
will be appreciated, the stenosis to be detected may 
range anywhere from a minor reduction in an artery's so 
diameter to a complete blockage or occlusion. The 
ECG/BSM system 10 is designed to detect the pres- 
ence and magnitude of stenoses, as well as the artery 
or arteries involved. 

Stenosis is typically related to ischemia. In that 55 
regard, in the presence of stenosis, the reduction in 
blood flow, may eventually result in ischemia of the arte- 
riole and cardiac tissue. The ECG/BSM system 1 0 does 
not respond to the stenosis directly, but rather the effect 



that the attendant ischemia has on action potentials. 
Given the spatial sensitivity of the mapping process 
employed, it is possible for the system 10 to achieve the 
desired sensitivity. 

In that regard, the 12-lead transformation and 
QRST area feature extraction described above have 
been found suitable for detecting not only the existence 
of stenosis, but also for discriminating between the loca- 
tion of the stenosis at such locations as the right coro- 
nary artery (RCA), left circumflex artery (LCX), and left 
anterior descending artery (LAD). By way of illustration, 
the performance of an ECG/BSM system 10 in this 
application will be briefly reviewed. 

The clinically based statistical analysis used to pro- 
gram the system 10 to identify the appropriate coeffi- 
cients for use was based upon a clinical population 
including 200 patients with suspected coronary artery 
stenosis. Coronary angiography indicated that 184 of 
the patients exhibited more than a 50% narrowing of at 
least one artery (1 13 LAD, 88 LCX, and 109 RCA). Out 
of this population of 200, 92 patients exhibited single 
artery stenosis or disease, 58 had double artery dis- 
ease, 34 had triple vessel disease, 51 patients had 
experienced a previous myocardial infarction, and 16 
had normal coronaries. An age-matched normal popu- 
lation of 137 patients was also used. 

Data was collected from the population using a 
ECG/BSM system 10 employing ten electrodes and 
using the first twelve basis functions computed for the 
192-lead transformation of the twelve-lead QRST inte- 
gral data. The statistical package used to identify the 
particular coefficients suitable for evaluating the degree 
of stenosis in the various arteries of interest was the 
BMDP package identified above. As will be appreciated, 
this software was used to separately analyze the clinical 
data for each locality of stenosis. 

Of the twelve resultant coefficients, the clinically 
based statistical analysis indicates that the first, second, 
fourth, and twelfth coefficients are appropriate for evalu- 
ating RCA stenosis, the first, second, third, and twelfth 
coefficients are appropriate for evaluating LCX stenosis, 
and the first, second, third, ninth, eleventh, and twelfth 
coefficients are appropriate for evaluating LAD stenosis. 
As will be appreciated, it may also be possible to identify 
certain coefficients that are more successful in identify- 
ing larger as opposed to smaller stenoses. 

The ability of the ECG/BSM system 10 to locally 
detect and evaluate the different stenoses of the clinical 
population was then evaluated, with the following 
results: 
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Artery 


Sensitivity 


Specificity 


Predictive 
Accuracy 


LAD 


85% 


80% 


82% 


LCX 


78% 


84% 


82% 


RCA 


82% 


82% 


82% 


CAD 


82% 


81% 


81% 



The CAD value listed is a cumulative estimation of the 
existence of stenosis in any vessel, obtained, for exam- 
ple, not through the development of separate CAD coef- 
ficients, but rather through a combined weighting of the 
LCX and LAD coefficients. The advantage of this com- 
bined weighting is that it frequently allows a more accu- 
rate overall detection of stenosis to be achieved. This 
approach has been found effective at detecting and 
localizing stenoses as small as 20% to 40% of vessel 
diameter. 

The areas of the ROC curves generated for the dif- 
ferent localities evaluated by the system 10 are 0.871 
for LAD, 0.894 for LCX, 0.866 for RCA, and 0.896 for 
CAD. The performance available from exercise treadmill 
testing (ETT) which is the currently accepted clinical 
standard, is a sensitivity of 66% and a specificity of 
84%. Thus, the performance of the ECG/BSM system 
1 0 represents a significant enhancement over ETT. The 
relative simplicity of the system 1 0 interfaces, its nonin- 
vasive nature, and the elimination of exercise as a test 
condition are among other advantages the ECG/BSM 
system 1 0 has over ETT. 

One enhancement of the preceding CAD analysis 
involves the way in which the magnitude of stenosis can 
be evaluated. In that regard, thallium testing can be 
used to gather clinical data concerning stenosis over 
five different cardiac regions: lateral, inferior, apical, 
anterior, and septal. As will be appreciated, the best 
coefficients for evaluating disease in these regions can 
then be determined in accordance with the analysis 
depicted in FIGURE 8. 

In one arrangement, a canonical variable V is com- 
puted for each of these five regions and compared 
against a cutpoint as part of a normal versus diseased 
diagnosis for each region. The magnitude of the steno- 
sis present in the LAD, LCX, or RCA can then be judged 
depending upon the number of the more limited regions 
in which stenosis is present. In that regard, the anterior 
and lateral subregions are combined in the evaluation of 
the magnitude of LAD stenosis. The lateral and apical 
regions are used in the LCX evaluation, while the infe- 
rior and septal regions are used in the RCA evaluation. 
With two subregions involved in the analysis of each 
artery, the resultant magnitude output can be simply 
GREATER THAN 50%, if both subregions are diseased, 
and LESS THAN 50%, if only one subregion is dis- 
eased. 
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Another cardiac characteristic that appears to be 
particularly suitable for analysis with 12-1 ead trans- 
formed BSM data is EP inducibility. Inducibility basically 
refers to the ability of an electrical shock, applied to the 

5 heart of a patient who has had a history of cardiac 
arrhythmia, such as fibrillation, to induce that arrhyth- 
mia. Drugs can be administered to the patient to reduce 
his or her susceptibility to the arrhythmia and a meas- 
urement of the patient's inducibility then becomes a 

w measure of the drugs' efficacy. 

The evaluation of inducibility is performed in largely 
the same manner as described for stenosis detection 
above. Of course, the population used to construct the 
analytical software relied upon by the system 10 differs. 

75 The coefficients selected for further evaluation are one, 
three, five and eight. Otherwise, the twelve-lead trans- 
formation of QRST area data may still be used and 
processed as described previously. 

Another particularly suitable application for the 

20 ECG/BSM system 10 involves the distinguishing of 
reversibly ischemic tissue from necrotic or dead tissue. 
In that regard, conventional evaluations involve the 
injection of thallium into a patient's blood stream. The 
patient's body is then imaged immediately after injection 

25 using gamma radiation, as well as some interval of time 
later. Healthy tissue will be immediately perfused and, 
as a result, exhibit thallium presence in the initial image. 
Reversibly ischemic tissue will not have thallium present 
in the initial image, but some perfusion will occur prior to 

30 the later imaging. Necrotic tissue will exhibit no per- 
fusion either time. 

Unfortunately, the images generated are extremely 
difficult to evaluate by all but the most experienced prac- 
titioners. Further, the invasive nature of the thallium 

35 injection and gamma radiation limits the desirability of 
the evaluation and the equipment is extremely expen- 
sive and not widely available. By constructing the 
ECG/BSM system 10 to collect cardiac data and select 
and evaluate coefficients based upon a thallium-based 

40 clinical population, however, a noninvasive diagnosis 
can be quickly and simply made. 

The design of system 10 parallels that described 
above. Because the thallium testing offers relatively 
detailed regional information, however, the regional 

45 sensitivity or localizing ability of the system 10 can be 
enhanced. As discussed above in connection with the 
evaluation of stenosis, additional regional information 
may offer some insight into the quantification of disease. 
Although not discussed in detail, the ECG/BSM 

so system 10 could also be used to evaluate other cardiac 
anomalies. For example, the system 10 might be used 
to evaluate the risk of arrythmia due to abnormal distri- 
bution of recovery properties. Similarly, the system 10 
might be used to evaluate reperfusion or restenosis 

55 after thrombolytic therapy of percutaneous angioplasty. 
As will be appreciated, the operation of system 10 
can be varied in a multitude of ways. For example, the 
data collected initially can be selected with the use of 
alternative waveform templates and tolerances to allow 
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cardiac cycles containing spurious features of interest 
to be evaluated. Similarly, the data collected can be 
reviewed over a number of cardiac cycles and a given 
signal averaged during any interval or specific time dur- 
ing the processing or the various signals averaged 5 
together. In addition, as suggested, above, the portion 
of the waveform that is analyzed over each cycle can be 
altered as desired. Further, the manner in which the 
data may be compressed, expanded, and statistically 
interpreted can vary, along with the cardiac features to w 
be detected and analyzed. 

Regardless of its form, the ECG/BSM system 10 
offers a number of advantages over existing systems. In 
that regard, most currently accepted evaluation tools 
involve an invasive examination that may significantly 15 
limit the ability or desirability of routinely examining 
larger populations to diagnose cardiac conditions at an 
early time, when they are often most readily treated. 

Further, in comparison to existing electrocardio- 
graphic systems, the system 1 0 advantageously offers 20 
the spatial sensitivity available only with the emerging 
BSM systems including a substantially larger array of 
electrodes. As a result, the system again has a more 
widespread applicability. Further, the feature extraction 
techniques employed can be used with any mapping- 25 
type electrode layout to produce an enhanced evalua- 
tion of the electrical activity of the patient's heart. 

Those skilled in the art will recognize that the 
embodiments disclosed herein are exemplary in nature 
and that various changes can be made therein without 30 
departing from the scope of the invention. In that regard, 
and as was previously mentioned, the hardware of the 
system may be included in a single device or as a plu- 
rality of discrete components or systems. Further, it will 
be recognized that although nine leads of data are 35 
described in connection with the transformation process 
disclosed, the advantages of the invention can be 
achieved with some variation of the number of leads 
employed, below the number previously deemed suita- 
ble. Because of the above and numerous other varia- 40 
tions and modifications that will occur to those skilled in 
the art, the following claims should define the invention. 

Claims 

45 

1 . A system for determining whether a patient has a 
coronary disease, comprising: 

a first set of electrodes (L, R, F, G, p1 , p2, p3, 
p4, p5, p6) for receiving a first set of cardiac so 
signals from a patient, at less than twenty elec- 
trodes; 

a data conversion interface (30) connected to 
receive the first set of cardiac signals and to 55 
convert the first set of cardiac signals into a dig- 
ital form; 

memory means (34) for storing a function that 



computes a canonical variable using three or 
more statistically determined coefficients taken 
from a set of coefficients; 

a central processing unit (32) that is pro- 
grammed to: 

a) transform the digitized set of cardiac sig- 
nals into a greater number of digitized val- 
ues, said set of values representing a set 
of body surface mapping signals such that 
the first set of cardiac signals appear as if 
they had been received from the patient 
using a set of body surface mapping elec- 
trodes, at more than twenty electrodes, 
and said body surface mapping electrodes 
had received the whole set of body surface 
mapping signals; said transformation using 
a mathematical expectation operator and 
stored referential data; 

b) to expand the transformed set of values 
representing the body surface mapping 
signals as a linear combination of prede- 
termined spatial basis functions to deter- 
mine said set of coefficients as being the 
coefficients of the linear expansion; 

c) to compute a value for the canonical var- 
iable using the stored function; and 

d) to compare the value of the canonical 
variable to a cut point above which the 
canonical variable indicates that the 
patient has the coronary disease and 
below which the canonical variable indi- 
cates the patient does not have the dis- 
ease; and 

a display (38) that indicates whether the patient 
has the coronary disease, the display being 
activated by the central processing unit if the 
value of the value of the canonical variable is 
greater than the cut point. 

2. The system of claim 1 , wherein the first set of elec- 
trodes comprises a four-electrode limb set 
(L.R.F.G) and a six-electrode precordial set (p1 , p2, 
p3, p4, p5, p6) 

3. The system of claim 1 , wherein the set of body sur- 
face mapping electrodes comprises a set of 192 
body surface mapping electrodes. 

4. The system of claim 1 , wherein the first set of car- 
diac signals obtained from a patient comprise inte- 
grals of a patient's QRST waveform. 

5. The system of claim 1 , wherein the set of body sur- 
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face mapping signals is represented as a vector, 
and wherein the central processing unit (32) is fur- 
ther programmed to: 

e) normalise the vector to reduce the influence 5 
of non-cardiac factors on the magnitude of the 
vector by scaling the vector to a unit length. 

6. The system of claim 1 , wherein the basis functions 
are identified using a Karhunen-Loeve (KL) expan- w 
sion. 

7. The system of claim 1, wherein the function that 
computes the canonical variable using the three or 
more statistically determined coefficients taken 75 
from the set of coefficients is derived by calculating 
sets of coefficients for a group of patients having 
the coronary disease and sets of coefficients of a 
group of patients not having the coronary disease, 
and performing a statistical analysis on the sets of 20 
coefficients obtained from both groups to select 
which three or more coefficients from the sets of 
coefficients best predict the presence of the coro- 
nary disease. 

25 

8. The system of claim 7, wherein the statistical anal- 
ysis used is a linear analysis. 

9. The system of claim 7, wherein the statistical anal- 
ysis used is a logarithmic analysis. 30 

10. The system of claim 1, wherein the coronary dis- 
ease is stenosis of the right coronary artery, left cir- 
cumflex artery, and left anterior descending artery; 
necrosis, accessory pathway blockages, bundle 35 
blockages, ischema, and cardiac nerve disorders. 

1 1 . The system of claim 1 , wherein the predetermined 
set of basis functions are the eigenvectors of the 
covariance matrix of the set of body surface map- 40 
ping signals and the three or more statistically 
determined coefficients used to compute the 
canonical variable are selected from the first twelve 
coefficients obtained by multiplying the set of body 
surface mapping signals by the predetermined set 45 
of basis functions. 

PatentansprOche 

1. System zur Feststellung, ob ein Patient an einer so 
Herzkranzgefa3krankheit leidet, umfassend: 

einen ersten Satz von Elektroden (L, R, F, G, 
p1, p2, p3, p4, p5, p6) zum Empfangen eines 
ersten Satzes von Herzsignalen von einem 55 2. 
Patienten an weniger als zwanzig Elektroden; 

eine Datenwandlerschnittstelle (30). die so ver- 
schaltet ist, da 3 sie den ersten Satz von Herz- 



signalen empfangt und den ersten Satz von 
Herzsignalen in digitate Form umwandelt; 

eine Speichereinrichtung (34) zum Speichern 
einer Funktion, welche unter Verwendung von 
drei Oder mehr statistisch bestimmten Koeffizi- 
enten aus einem Satz von Koeffizienten eine 
kanonische Variable berechnet; 

eine zentrale Verarbertungseinheit (32), die 
programmiert ist, urn 

a) den digitalen Satz von Herzsignalen 
umzuwandeln in eine grOBere Anzahl digi- 
talisierter Werte, wobei der Satz von Wer- 
ten einen Satz von KGrperoberflachen- 
Kartografiesignalen reprasentiert. derart, 
da 13 der erste Satz von Herzsignalen so 
erscheint, als sei er von dem Patienten 
unter Verwendung eines Satzes von KOr- 
peroberflachen-Kartografieelektroden an 
mehr als zwanzig Elektroden aufgenom- 
men worden, und die KCrperoberflachen- 
Kartografieelektroden den gesamten Satz 
von Korperoberflachen-Kartografiesigna- 
len empfangen hatten; wobei die Umwand- 
lung einen mathematischen Erwartungs- 
operator und gespeicherte Referenzdaten 
verwendet; 

b) den umgewandelten Satz von Werten, 
der die KOrperoberflachen-Kartografiesi- 
gnale als Linearkombination vorbestimm- 
ter raumlicher Basisfunktionen reprasen- 
tiert, zu erweitern, urn den Satz von Koeffi- 
zienten zu bestimmen als Koeffizienten 
der linearen Erwerterung, 

c) einen Wert fur die kanonische Variable 
unter Verwendung der gespeicherten 
Funktion zu berechnen; und 

d) den Wert der kanonischen Variablen zu 
vergleichen mit einem Schnittpunkt, ober- 
halb dessen die kanonische Variable 
bedeutet, daf3 der Patient die Herzkranz- 
gefaGkrankheit aufweist, und unterhalb 
dessen die kanonische Variable angibt, 
da 3 der Patient nicht an der Krankheit lei- 
det; und 

eine Anzeige (38), die angibt, ob der Patient an 
der HerzkranzgefaBkrankheit leidet, wobei die 
Anzeige durch die zentrale Verarbeitungsein- 
heit aktiviert wird, falls der Wert der kanoni- 
schen Variablen grOBer ist als der Schnittpunkt. 

System nach Anspruch 1 , bei dem der erste Satz 
von Elektroden einen vier Elektroden umfassenden 
Extremitatensatz (L, R, F, G) und einen sechs Elek- 
troden umfassenden Prakordialsatz (p1 , p2, p3, p4, 
p5, p6) umfaBt. 
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3. System nach Anspruch 1, bei dem der Satz von 
KOrperoberflachen-Kartografieelektroden einen 
Satz von 192 Kdrperoberflachen-Kartografieelek- 
troden umfaBt. 

5 

4. System nach Anspruch 1 , bei dem der erste Satz 
von Herzsignalen, die von einem Patienten erhalten 
werden, Integrale einer QRST-Wellenform des Pati- 
enten umfaBt. 

10 

5. System nach Anspruch 1, bei dem der Satz von 
Korperoberf lachen-Kartografiesignalen als ein Vek- 
tor dargestelH wird, wobei die zentrale Verarbei- 
tungseinheit (32) auGerdem programmiert ist, urn: 

15 

e) den Vektor zu normieren und so den EinfluB 
von herzfremden Faktoren auf den Betrag des 
Vektors zu verringern, indem der Vektor auf 
eine Einheitslange skaliert wird. 

20 

6. System nach Anspruch 1 , bei dem die Basisfunktio- 
nen unter Verwendung einer Karhunen-Loeve-(KL- 
JErweiterung identifiziert werden. 

7. System nach Anspruch 1, bei dem die Funktion, 25 
welche die kanonische Variable unter Verwendung 
von drei Oder mehr statistisch bestimmten Koeffizi- 
enten aus dem Satz von Koeffizienten berechnet, 
dadurch abgeleitet wird, daB Satze von Koeffizien- 
ten for eine Gruppe von Patienten, welche die Herz- 30 
kranzgefaBkrankheit aufweisen, und Satze von 
Koeffizienten einer Gruppe von Patienten, die die 
Herz kranzgefaBkrankheit nicht aufweisen, berech- 
net werden, und eine statistisch e Analyse der 
Satze von Koeffizienten von beiden Gruppen 35 
durchgefuhrt wird, urn herauszufinden, welche drei 
Oder mehr Koeffizienten aus den Satzen von Koef- 
fizienten das Vorliegen der HerzkranzgefaBkrank- 
heit am besten vorhersagen. 

40 

8. System nach Anspruch 7, bei dem die verwendete 
statistisch e Analyse eine lineare Analyse ist. 

9. System nach Anspruch 7, bei dem die verwendete 
statistische Analyse eine logarithmische Analyse 45 
ist. 

10. System nach Anspruch 1, bei dem die Herzkranz- 
gefaBkrankheit eine Verengung der rechten Koro- 
nararterie, der linken Zircumflexarterie und der so 
linken vorderen Abwartsarterie; Nekrose, Neben- 
wegblockierungen, Bundelblockagen, Ischamie 
und HerznervenstOrungen umfaBt. 

11. System nach Anspruch 1, bei dem der vorbe- 55 
stimmte Satz von Basisfunktionen die Eigenvekto- 
ren der Kovarianzmatrix des Satzes von 
Korperoberflachen-Kartografiesignalen umfaBt, 
und die drei Oder mehr statistisch bestimmten Koef- 



fizienten zum Berechnen der kanonischen Varia- 
blen ausgewahlt sind aus den ersten zwGlf 
Koeffizienten, die erhalten werden durch Muttipli- 
zieren des Satzes von KOrperoberflachen-Karto- 
grafiesignalen mit dem vorbestimmten Satz von 
Basisfunktionen. 

Revendications 

1. Systeme pour determiner si un patient souffre une 
maladie coronarienne, comprenant : 

un premier ensemble d'electrodes (L,R,F,G,- 
p1,p2, p3.p4,p5,p6) pour recevoir un premier 
ensemble de signaux cardiaques emanant d'un 
patient, au niveau de moins de vingt electro- 
des; 

une interface de conversion de donnees (30) 
connectee de maniere a recevoir le premier 
ensemble de signaux cardiaques et a convertir 
le premier ensemble de signaux cardiaques 
sous une forme numerique; 
des moyens de memoire (34) pour memoriser 
une fonction qui calcule une variable canoni- 
que utilisant trois ou un plus grand nombre de 
coefficients determines statistiquement, tires 
d'un ensemble de coefficients; 
une unite centrale de traitement (32), qui est 
programmed pour : 

a) transformer I'ensemble numerise de 
signaux cardiaques en un nombre plus 
eleve de valeurs numerisees, ledit ensem- 
ble de valeurs representant I'ensemble de 
signaux d'etablissement d'une carte de la 
surface du corps de sorte que le premier 
ensemble de signaux cardiaques apparatt 
comme si ces signaux avaient ete recus 
par le patient en utilisant un ensemble 
d'electrodes d'etablissement de la carte de 
la surface du corps, au niveau de plus de 
vingt electrodes, et comme si lesdites elec- 
trodes d'etablissement de la carte de la 
surface du corps avaient recu I'ensemble 
complet de signaux d'etablissement de la 
carte de la surface du corps; ladite trans- 
formation utilisant un operateur d'espe- 
rance mathematique, et des donnees de 
reference memorisees; 

b) developper rensemble tran storm e de 
valeurs representant les signaux d'etablis- 
sement de la carte de la surface du corps, 
sous la forme d'une combinaison lineaire 
de fonctions spatial es de base predetermi- 
nes pour determiner ledit ensemble de 
coefficients comme etant les coefficients 
de developpement lineaire; 

c) calculer une valeur pour la variable 
% canonique en utilisant la fonction memori- 
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see; et 

d) comparer la valeur de la variable cano- 
nique a un point de separation, au-dessus 
duquel la variable canonique indique que 
le patient possede la maladie corona- 5 
rienne et au-dessous duquel la variable 
canonique indique que le patient ne pre- 
sent e pas cette maladie; et 

un dispositrf d'affichage (38) qui indique si le w 
patient possede la maladie coronarienne, raff i- 
chage etant active par I'unite centrale de traite- 
ment si la valeur de la variable canonique est 
superieure au point de separation. 

15 

Systeme selon la revendication 1, dans lequel le 
premier ensemble d'6Iectrodes comprend un 
ensemble de branches a quatre electrodes 
(L.R.F.G) et un ensemble precordial de six eiectro- 
des (p1 ,p2,p3,p4,p5,p6). 20 

Systeme selon la revendication 1, dans lequel 
Tensemble d'electrodes reproduisant I'image de la 
surface du corps comprend un ensemble de 192 
Electrodes d*6tablissement de la carte de la surface 25 
du corps. 

Systeme selon la revendication 1, dans lequel le 
premier ensemble de signaux cardiaques obtenus 
chez un patient comprend les integrates d'une 30 
forme d'onde QRST du patient. 



quels trois coefficients ou plus tires des ensembles 
de coefficients deer i vent le mieux la presence de la 
maladie coronarienne. 

8. Systeme selon la revendication 7, dans lequel 
I'analyse statist ique utilis&e est une analyse 
lineaire. 

9. Systeme selon la revendication 7, dans lequel 
I'analyse statistique utilisee est une analyse loga- 
rithmique. 

10. Systeme selon la revendication 1, dans lequel la 
maladie coronarienne est une stenose de I'artere 
coronaire droite, de I'artere circonf lexe gauche et 
de I'artere descendante anterieure gauche; une 
necrose, des blocages de voies accessoires, des 
blocages de faisceaux, Tischemie et des troubles 
cardiaques nerveux. 

11. Systeme selon la revendication 1, dans lequel 
I'ensemble predetermine de fonctions de base sont 
les vecteurs propres de la matrice de covariance de 
I'ensemble des signaux reproduisant I'image de la 
surface du corps, et les trois coefficients ou plus, 
determines statistiquement, utilises pour calculer la 
variable canonique sont selection nes a parti r des 
douze premiers coefficients obtenus en multipliant 
i'ensemble de signaux d'etablissement de la carte 
de la surface du corps par I'ensemble predetermine 
de fonctions de base. 



Systeme selon la revendication 1, dans lequel 
I'ensemble de signaux reproduisant I'image de la 
surface du corps est represents par un vecteur, et 35 
I'unite centrale de trartement (32) est en outre pro- 
grammee pour : 



e) normaliser le vecteur pour recJuire Pinf luence 
de facteurs non cardiaques sur I'amplitude du 40 
vecteur par reglage du module du vecteur sur 
une longueur unite. 

Systeme selon la revendication 1 , dans lequel les 
fonctions de base sont identrfi6es moyennant I'utili- 45 
sation d'un d6veloppemerrt de Karhunen-Loeve 
(KL). 



Systeme selon la revendication 1, dans lequel la 
fonction, qui calcule la variable canonique en utili- so 
sant trois ou un plus grand nombre de coefficients 
determines statistiquement et tires de I'ensemble 
de coefficients, est obtenue par calcul d'ensembfes 
de coefficients pour un groupe de patients pr6sen- 
tant ta maladie coronarienne et d'ensembles de 55 
coefficients d'un groupe de patients ne presentant 
pas la maladie coronarienne, et en executant une 
analyse statistique des ensembles de coefficients 
obtenus a partir de groupes pour seiectionner les- 
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